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Abstract. We present first measurements of charged and neutral particle-flow correlations in pp colli- 
sions using the ATLAS calorimeters. Data were collected in 2009 and 20f0 at centre-of-mass energies of 
900 GeV and 7 TeV . Events were selected using a minimum- bias trigger which required a charged particle 
in scintillation counters on either side of the interaction point. Particle flows, sensitive to the underlying 
event, are measured using clusters of energy in the ATLAS calorimeters, taking advantage of their fine 
granularity. No Monte Carlo generator used in this analysis can accurately describe the measurements. 
The results are independent of those based on charged particles measured by the ATLAS tracking systems 
and can be used to constrain the parameters of Monte Carlo generators. 

PACS. f 2.38.Qk Experimental tests - 12. 38. Aw General properties of QCD 



1 Introduction 

All hard parton-parton interactions in pp collisions are ac- 
companied by additional processes which collectively con- 
tribute additional particles to those from the hard scat- 
ter and which are termed the underlying event (UE). It 
is impossible to uniquely separate the UE from the hard 
scattering process on an event-by-event basis. However, 
observables can be measured which are sensitive to prop- 
erties of the UE. In order to make high-precision measure- 
ments, the UE must be modelled using phenomenological 
models in Monte Carlo generators pQ. Such models must 
be tuned to experimental data. In the past, such studies 
have only been performed using tracks [2JJ5]. 

Many physics processes to be studied with the AT- 
LAS detector [6] require precision measurements of jets 
and missing transverse energy obtained principally from 
the calorimeter system. Therefore, it is important that the 
UE measurements are performed using the same instru- 
mental environment and reconstructed objects as those 
for the calorimeter-based measurements. The fine granu- 
larity of the ATLAS calorimeter allows the definition of 
three-dimensional clusters of energy which are closely as- 
sociated with individual particles 018] . 

A study of the UE using charged-track densities was 
recently performed by ATLAS [5]. The present paper ex- 
tends this measurement by reconstructing particle densi- 
ties using calorimeter clusters in the region which is most 
sensitive to the soft QCD processes responsible for the 
UE; the "transverse" region as shown in Fig. [1] The az- 
imuthal angular distance between a leading particle in 
transverse momentum (pt) and other particles is given by 



leading particle 




Fig. 1. A schematic representation of regions in the azimuthal 
angle <j> with respect to the leading particle (shown with the 
arrow). In this analysis, the leading particle corresponds to the 
cluster with the largest transverse momentum. 



A(f> = <fi — 0ioad, where <f> is the azimuthal angle of a parti- 
cle and (^i C ad is the azimuthal angle of the leading particle. 
The transverse region, defined as 60° < \A(j>\ < 120°, is 
most sensitive to the UE since it is perpendicular to the 
axis of hardest scattering, approximated by the direction 
of the leading particle. As is the case for charged particles, 
the number density of the clusters and their transverse en- 
ergy density in this region are sensitive, discriminating ob- 
servables for UE studies. These distributions are corrected 
for detector effects to give a measure of the particle ac- 
tivity in the UE and to provide new characteristics which 
can be used to tune models included in Monte Carlo gen- 
erators. 
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The analysis using calorimeter clusters has several im- 
portant features. Firstly, its results are sensitive to the en- 
tire hadronic final state, including neutral particles, which 
constitute about 40% of all produced particles. Secondly, 
the analysis based only on calorimeter clusters has com- 
pletely independent experimental uncertainties compared 
to the corresponding analysis [5] using charged particles. 
Finally, as discussed earlier, since jet reconstruction is 
based almost entirely on energy deposition in the calorime- 
ter, the results of this UE analysis can be used directly to 
estimate the effect of the underlying event on any jet- 
based measurement. 



2 The ATLAS Detector 

The ATLAS detector [6] at the Large Hadron Collider was 
designed to study a wide range of physics. It covers almost 
the entire solid angle around the collision point with layers 
of tracking detectors, calorimeters and muon chambers. 

Charged tracks and vertices are reconstructed with the 
inner detector which consists of a silicon pixel detector, a 
silicon strip detector and a transition radiation tracker, all 
immersed in a 2 tesla magnetic field provided by a super- 
conducting solenoid. For the measurements presented in 
this paper, the high-granularity calorimeter systems are 
of particular importance. The ATLAS calorimeter system 
provides fine-grained measurements of shower energy de- 
positions over a large range in pseudorapiditjQ- Electro- 
magnetic calorimetry in the range < 3.2 is provided by 
liquid argon (LAr) sampling calorimeters. This calorime- 
ter system provides measurements of the shower energy 
in up to four depth segments and with transverse gran- 
ularity that ranges from 0.003 x 0.10 to 0.05 x 0.025 in 
Sr] x 8<j), depending on depth segment and rapidity. The 
hadronic calorimetry in the range |7y| < 1.7 is provided by 
a steel/scintillator-tile sampling calorimeter. This system 
provides measurements of the shower energy deposition in 
three depth segments at a transverse granularity of typi- 
cally 0.1 x 0.1. In the end-caps (\r)\ > 1.5), LAr technol- 
ogy is used for the hadronic calorimeters that match the 
outer rj limits of the end-cap electromagnetic calorimeters. 
This system provides four measurements in depth of the 
shower energy deposition at a transverse granularity of ei- 
ther 0.1x0.1 (1.5 < |ry| < 2.5) or0.2x0.2 (2.5 < < 3.2). 
The LAr forward calorimeters provide both electromag- 
netic and hadronic energy measurements and extend the 
calorimeter coverage from = 3.2 to \r)\ = 4.9. To mea- 
sure the energy of photons and electrons, all calorimeter 

The ATLAS reference system is a Cartesian right-handed 
co-ordinate system, with the nominal collision point at the ori- 
gin. The anti-clockwise beam direction defines the positive z- 
axis, while the positive a;-axis is defined as pointing from the 
collision point to the centre of the LHC ring and the posi- 
tive y-axis points upwards. The azimuthal angle <f> is measured 
around the beam axis, and the polar angle 8 is the angle mea- 
sured with respect to the z-axis. The pseudorapidity is given 
by r] = — lntan(0/2). Transverse momentum is defined relative 
to the beam axis. 



cells are initially calibrated to the electromagnetic energy 
scale using test-beam data [9Tfl4]. 

This analysis is based on the properties of topological 
clusters in the calorimeter, which represent an attempt 
to reconstruct three-dimensional energy depositions as- 
sociated with individual particles 015] ■ The topological- 
cluster algorithm proceeds through the following steps. 
Nearest neighbours are collected around seed cells, which 
are cells with an absolute signal greater than 4er above the 
noise level [SHE]- Then, neighbouring cells are collected 
into the cluster if the absolute value of their signal sig- 
nificance is above a secondary seed threshold of 2a. All 
surrounding cells are iteratively added to the cluster until 
no further secondary seeds are among the direct neigh- 
bours. A final analysis of the resulting cluster looks for 
multiple local signal maxima; in the case of more than 
one maximum in a given cluster, it is split into smaller 
clusters along the signal valleys between the maxima. 

The analysis presented in this paper uses calibrated 
topological clusters [8j. The clusters are classified as re- 
lated to electromagnetic or hadronic energy deposits, us- 
ing detailed information on the cluster topology. Then, 
the reconstructed cluster energy is corrected for the non- 
compensating nature of the ATLAS calorimeter and for 
inactive material. 



3 Data selection 

The data taken at y/s = 900 GeVwere collected during 
6-15 December 2009. During this running period, there 
were approximately 3% non-functional channels in the tile 
hadronic calorimeter and approximately 1% non- functional 
channels in the LAr calorimeters [121 El- For an inte- 
grated luminosity of 7 /ib _1 , a total of 455 thousand events 
were collected from colliding proton bunches in which a 
minimum-bias trigger recorded one or more hits in the 
scintillators on either side of the detector. 

The events to be analysed were selected using a pro- 
cedure identical to that described in Refs. [5HTS]. Events 
were required to have a primary vertex which is recon- 
structed using at least two tracks with transverse mo- 
menta pt > 100 MeV and a transverse distance of closest 
approach with respect to the beam-spot position |16j of 
less than 4 mm. 

This analysis uses topological clusters with > 0.5 
GeVand < 2.5 in order to have the same kinematic 
rangeQ as for the previous UE study based on tracks [S]. 
Additional selection criteria were applied to select good- 
quality clusters: 1) to reject the cosmic and noise back- 
ground, the leading cell energy of the cluster is required 
to be less than 90% of the cluster energy; 2) the energy 
sampling maximum must be in a calorimeter region with 
good calibration; 3) the fraction of energy in the cluster 
associated with problematic cells (or dead cells where the 

2 The topological clusters are treated as massless particles, 
and we choose to refer to both clusters and stable particles in 
terms of px- The same symbol pt is also used to represent the 
track transverse momentum. 
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energy contribution is obtained by interpolation from ad- 
jacent cells) should be less than 50%. 

Data at = 7 TeVwere collected between 30 March 
and 27 April 2010. Only a fraction of the 7 TeVdata, cor- 
responding to an integrated luminosity of about 230 ^b -1 , 
was used. In total, about 7.7 million events were analysed. 
Event selection was similar to that for the 900 GeV data, 
but included the additional requirement of a single pri- 
mary vertex [5"lll5] to remove events containing more than 
one pp interaction. 

4 Monte Carlo simulation 

The QCD predictions for the hadronic final state in in- 
elastic pp collisions are based on several Monte Carlo gen- 
erators. The PYTHIA 6.4 Monte Carlo generator [17] is 
used as the primary generator for comparisons with the 
data. The MC09 tune [H] of this model was performed 
by ATLAS. It uses the pr-ordered parton shower with the 
MRST LO* parton-density function [T5], followed by frag- 
mentation into final-state particles using the Lund string 
model [2U] . The parameters of this generator were adjusted 
to describe charged-particle multiplicity distributions in 
minimum-bias events measured at y/s — 630 GeV and 
y/s = 1.8 TeVin pp collisions [2"T] . Diffractive processes 
are not included in the simulation for the main samples, 
but were used for systematic checks (Sect. [7]). In addition 
to the MC09 tune, the following two PYTHIA parame- 
ter sets are also used: (1) the PerugiaO set [H] in which 
the soft-QCD part is tuned using only minimum-bias data 
from the Tevatron and CERN pp colliders; (2) the DW [23] 
PYTHIA tune, which uses virtuality-ordered showers and 
was derived to describe the CDF Run II underlying event 
and Drell-Yan data. 

The data are also compared to the PHOJET Monte 
Carlo generator [24] . which includes a simulation of the 
diffractive component. This generator is based on the two- 
component Dual Parton Model which includes soft hadronic 
processes described by Pomeron exchange and semi-hard 
processes described by perturbative parton scattering. The 
description of the fragmentation is the same as in the 
PYTHIA generator. 

In addition, the HERWIG Monte Carlo generator [25] 
126) was used for comparisons with the data. This gener- 
ator has similar matrix-element calculations as PYTHIA, 
but uses the cluster fragmentation model to hadronise par- 
tons into hadrons. HERWIG is interfaced with the JIMMY 
model [27] in order to describe multiple parton interac- 
tions. 

Monte Carlo events were processed through the AT- 
LAS detector simulation program [55], which is based on 
Geant4 [53]. They were reconstructed using the same 
trigger and event selection as for the data. The size and po- 
sition of the collision beam-spot and the detailed descrip- 
tion of detector conditions during the data-taking runs 
were included in the simulation. 

Monte Carlo events after the detector simulation pro- 
gram were used for correcting the data to the stable- 
particle level defined as follows. The PYTHIA MC09 is 



used to generate the primary samples for unfolding the 
effects of the detector. Monte Carlo stable particles are 
selected if their mean lifetimes are larger than 3 • 10~ n 
seconds. Neutrinos are excluded from consideration. Ac- 
cording to this definition, Kg, A and are among those 
treated as stable particles. This definition allows a direct 
comparison between the results of previous track-based 
studies [5] and the present measurement. 



5 Properties of calorimeter clusters 

Figures [2] and [3] show the distributions of pr and rj for 
topological clusters in data and simulated PYTHIA MC09 
events at y/s = 900 GeV and y/s = 7 TeV , respectively. 
The distributions in each case are normalised to the num- 
ber of entries. In addition, the ratio plots show the ratio 
of simulation to data in the transverse region alone. The 
figures show overall good agreement between the data and 
the PYTHIA MC09 tune, with 20% discrepancies in some 
phase-space regions. While not shown in these figures, the 
PerugiaO tune agrees with the data to a similar extent. 
The contribution of the discrepancy in the high-px tail 
is expected to be small on particle densities measured 
at pt > 0.5 GeV , and it was taken into account using 
re-weighting as described below. The observed differences 
between the data and the PYTHIA MC09 event sample 
for the r\ distributions are addressed in the studies of sys- 
tematic uncertainties. 

FiguresUJa) and[SJa) show the multiplicity of topolog- 
ical clusters, with pt > 0.5 GeV and \rj\ < 2.5, versus the 
number of stable particles (charged and neutral) in simu- 
lated events for yfs = 900 GeV and ^fs — 7 TeV . A strong 
correlation is observed between the number of topological 
clusters and the number of stable particles, indicating that 
clusters are a good representation of the particle activity 
in inelastic pp events. 

Figures SJb) and [5th) show the correlation between 
the number of topological clusters and the number of pri- 
mary tracks selected in the same way as in the track-based 
studies [5][TS]. These figures also show a strong correla- 
tion. The Monte Carlo simulation shown in Fig. HJc) and 
[5jc) reproduces these distributions well: the means and 
the root-mean-square deviations of one-dimensional pro- 
jections of these distributions agree with the Monte Carlo 
simulation within less than one percent for N (tracks) > 4. 
For events with a lower number of tracks, the data show a 
smaller mean value of the projection onto the a;- axis than 
seen in the PYTHIA MC09 simulation. This is attributed 
to the absence of diffraction in the generated samples. 

A Monte Carlo simulation study based on PYTHIA 
MC09 indicates that the probability that a second parti- 
cle lies within y/ 8rf + 8<p 2 < 0.2 of a first in the selected 
inelastic pp events is below 1%. This simplifies the present 
analysis since there is negligible potential bias due to clus- 
ter overlap. 

For the UE studies based on topological clusters, a 
good position measurement is required. The quality of the 
position reconstruction of the clusters was studied by com- 
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Fig. 2. A comparison between uncorrected data and the Monte Carlo simulation for topological cluster pt (a) and rj (b) for 
pp collisions at y/s = 900 GeV . The ratio plots show the inclusive sample (solid lines) and the transverse region (dashed lines). 
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Fig. 3. A comparison between uncorrected data and the Monte Carlo simulation for topological cluster px (a) and rj (b) for 
pp collisions at y/s — 7 TeV . The ratio plots show the inclusive sample (solid lines) and the transverse region (dashed lines). 



paring the impact point of charged particles with the as- 
sociated cluster position in the calorimeter. Charged par- 
ticles are deflected in the magnetic field of the solenoid. 
Their trajectories are extrapolated to the calorimeter us- 
ing a Monte Carlo simulation which includes a detailed 
field map as well as the effect of the material in front of 
the calorimeter. The Monte Carlo simulation describes the 
topological-cluster positions relative to the positions of the 
extrapolated tracks on the surface of the LAr calorimeter 



within the granularity of its second layer (0.025 x 0.025 in 
5rj x 54>). 

As the correction for detector effects is based on the 
Monte Carlo simulation, an essential issue is the accu- 
racy with which the simulation reproduces the energy re- 
construction in the calorimeter. For charged particles, the 
energy scale was studied [30.31 using isolated tracks by 
extrapolating tracks to the calorimeter surface and match- 
ing them to topological clusters. The average value of the 
ratio E/p was reconstructed, where E is the cluster energy 
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Fig. 4. ATLAS data at y/s = 900 GeV : The correlations between the multiplicities of, (a) topological clusters (TV(clusters)) and 
stable particles (AT (truth)) from simulated pp interactions, (b) topological clusters and primary reconstructed tracks (TV (tracks)) 
from pp interactions, and (c) topological clusters and primary reconstructed tracks from simulated pp interactions. Inelastic 
events generated by PYTHIA MC09 (without diffraction) passed through the selection were used to produce the plots (a) and 
(c). 




Fig. 5. ATLAS data at y/s = 7 TeV : The correlations between the multiplicities of, (a) topological clusters (TV(clusters)) and 
stable particles (TV (truth)) from simulated pp interactions, (b) topological clusters and primary reconstructed tracks (TV(tracks)) 
from pp interactions, and (c) topological clusters and primary reconstructed tracks from simulated pp interactions. Inelastic 
events generated by PYTHIA MC09 (without diffraction) passed through the selection were used to produce the plots (a) and 
(c). 



in the calorimeter and p is the track momentum. Figure [6] 
shows the average response (E/p) for calibrated topologi- 
cal clusters as a function of r] for tracks with p > 0.5 GeV . 
The data and PYTHIA MC09 agree within 5% in most f] 
regions, while discrepancies increase in the transition re- 
gion (1.5 < \r)\ < 1.8) between barrel and end-cap. 

To estimate the relative energy-scale uncertainty, the 
double ratio (E jp) mc I '(E '/p) was calculated, where the 
ratio (E/p)mc was determined from the Monte Carlo sim- 
ulation. The double ratio as a function of r/ is shown in 
Fig. H] (bottom). The double-ratio distributions were mea- 
sured for a wide range of track momenta and rj as de- 
scribed in Refs. [301111]. 



The comparison between data and Monte Carlo pre- 
dictions for the shapes of the E/p distribution is shown 
in Fig. The peak at zero corresponds to isolated tracks 
that have no associated cluster in the calorimeter. These 
are predominantly due to hadronic interactions in the ma- 
terial in front of the calorimeter [3T]. The contribution of 
the discrepancies observed for E/p — between the data 
and PYTHIA MC09 to uncertainties on the reconstruc- 
tion efficiencies of topological clusters is below 1%. This 
effect was taken into account as described in Sect. [7] More 
details on the energy scale of topological clusters can be 
found in Refs. [3j[^]. 

The energy scale for electromagnetic clusters was es- 
timated using the 7r° peak reconstructed in inelastic pp 
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Fig. 6. The average E/p in different r\ bins for isolated topolog- 
ical clusters matched to charged tracks in inelastic pp events at 
y/s — 7 TeV for track momentum p larger than 0.5 GeV . A sim- 
ilar level of agreement between data and the Monte Carlo sim- 
ulation was obtained for the y/s = 900 GeV data (not shown) . 



events. The selection criteria for calibrated topological 
clusters were the same as for the present analysis. The 7r 
peak positions for data and PYTHIA MC09 agree within 
3% for all T) regions. 

Correction of the observed distributions to the parti- 
cle level requires a reliable description of the cluster mul- 
tiplicity distribution by the simulated event sample. This 
was studied by examining cluster multiplicities in bins of 
track multiplicity using projections of the two-dimensional 
distributions shown in Figs. 0] and [SJ The observed differ- 
ences are propagated into the systematic uncertainties as 
discussed in Sect. [7J 

For the UE studies, the so-called "leading" clusters, 
i.e. clusters with the largest transverse momenta, pl^ , 
are used to define an event orientation. Such clusters are 
typically inside the most energetic jets. To verify this, jets 
were reconstructed with the anti-k t algorithm |33j with a 
distance parameter of 0.4, a minimum p^ requirement of 
5 GeV and \rj\ < 2.5. Then, the distance in rj-<p between the 
leading topological cluster and the centre of the leading jet 
was calculated. It was shown that, in the vast majority of 
cases, the leading cluster is inside a leading jet with only a 
small fraction (~ 10%) of clusters opposite the leading jet 
in (j). This feature is well reproduced by the Monte Carlo 
simulation. 



To verify the Monte Carlo performance for , the 
ratio of pl£ ad of topological clusters to pl£ ad of primary 
tracks was reconstructed. The agreement between the data 
and the PYTHIA MC09 tune for such distributions was 
found to be within ±5% in most regions, while discrep- 
ancies at the level of 20% were found for the tails of the 
ratio distributions. The impact of such discrepancies in the 
simulation of the pl£ ad resolution on the final measurement 
has been estimated as discussed in Sect. [7] 

Monte Carlo studies show that the rate of events in 
which a low-px particle is reconstructed as a high-px clus- 
ter is not negligible. This results in a low purity for topo- 
logical clusters at high pl^ ad . Therefore, the analysis was 
performed for leading topological clusters with transverse 
momenta less than 8 GeV (14 GeV ) for the = 900 GeV 
(7 TeV ) data in order to limit this effect and to ensure that 
the reconstruction purity even at the highest transverse 
momenta considered is larger than 50%. 

6 Measured observables and correction 
procedure 

Following earlier track- based analyses [5], particle densi- 
ties are studied as a function of the distance A<j> in the 
azimuthal angle between the leading cluster and all other 
clusters in an event, and as a function of pt of the leading 
cluster in the event. The scalar pt sum for stable particles 
per unit area in rj— <j> in the transverse region is also pre- 
sented. This provides complementary information to that 
which can be obtained from the particle densities. 

The particles and clusters are required to have pt > 
0.5 GeV and \rj\ < 2.5. Clusters are selected if they pass 
the criteria described in Scct.[3J The measured observables 
at the particle and detector levels are: 

— pl^ ad - Transverse momentum of the stable particle 
with maximum m the event. At the detector level, 
this corresponds to the transverse momentum of the 
selected topological cluster with maximum p^ in the 
event. 

— d(N) I dA(f> - The average number of particles as a func- 
tion of the azimuthal-angle difference between the lead- 
ing particle and other particles in an event. The lead- 
ing particle at A<p = is excluded from this distri- 
bution. At the detector level, it corresponds to the 
mean number of selected topological clusters as a func- 
tion of the azimuthal-angle distance between the lead- 
ing topological cluster and other clusters in an event. 
This density [5] is defined per unit of pseudorapidity 
as N/(N ev ■ (?7 max — ??min)), where N is the number of 
entries in Acj) bins, i] max — i] mm — 5 represents the full 
pseudorapidity range, and N ev is the number of events 
selected by requiring a particle with pl^ ad above the 
specified value. 

— (d 2 N / drjdcj)) - Mean number of stable particles per unit 
area in r)-<p. At the detector level, this corresponds to 
the mean number of selected topological clusters per 
unit area in 77— 4>. This density is measured as a function 
of p l f ad [5]. 
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Fig. 7. .E/p distributions at -^/s = 7 TeV for topological clusters matched to tracks in several bins of track momentum: (a) 
0.5 < p < 1.2 GeV , (b) 1.2 < p < 2.2 GeVand (c) 2.2 < p < 10 GeV . The peak at zero corresponds to the events without a 
good match between a topological cluster and a track. A similar level of agreement between data and Monte Carlo was obtained 
for the y/s = 900 GeV data (not shown). 



— (d 2 52 Pt/ drjdcf)) - Mean scalar px sum for stable par- 
ticles per unit area in rj—cf>. At the detector level, this 
corresponds to the mean scalar px sum for selected 
topological clusters per unit area in r)-(f>. This quantity 
is defined following the convention used in the previous 
ATLAS publication [5]. 

A bin-by-bin correction procedure is used to correct 
the observed distributions to the stable-particle level. The 
correction factors 

A scn 

^ = ^dct ' 

are evaluated separately for each observable. In the above 
expression, A gcn is calculated at the stable-particle level 
of PYTHIA MC09 and .4 dot is calculated after full detec- 
tor simulation and reconstruction. The corrected value for 
an observable is found by multiplying its measured value 
by the relevant correction factor C. These factors correct 
the data to the stable-particle level and include the effects 
of event selection, reconstruction efficiency, bin migrations 
and smearing, including the case when the leading parti- 
cle is mis-identified and a cluster corresponding to a sub- 
leading particle is used to define the event orientation and 

„lcad 
Pt ■ 

The bin-by-bin correction depends on the choice of the 
Monte Carlo event generator. This affects the efficiency 
correction (mainly due to variations in particle types) and 
the purity (different stable-particle level distributions have 
different fractions of poorly reconstructed objects in each 
bin as well as different bin migrations). To reduce the 
model dependence of the correction procedure, bin-by- 
bin migrations were minimised by using bin sizes larger 
than the reconstruction resolutions for the distributions 



presented. In addition, the analysis was restricted to the 
p!j; ad ranges where the purity of leading clusters is above 
50% (see Sect. EJ). 

The bin- by-bin correction factors for the particle densi- 
ties typically have values of around 1.3 and do not exceed 
1.4. The largest single contributor is the reconstruction 
inefficiency of topological clusters, which leads to a bin- 
by-bin correction factor of approximately 1.2 on average 
and has a maximum value of 1.3 at low px- The other 
significant contributor is the event reorientation which re- 
sults from inefficiency of the reconstruction of the leading 
topological cluster in an event. This causes bin migrations, 
which were studied by replacing the leading cluster pJg ad 
by the leading track p!j5 ad (track), for which the efficiency 
is known to be high [5] . The bin- by-bin corrections for the 
average scalar px sum have a maximum value of 1.5 for 
low p^ ad and decrease to 1.3 for p!j5 ad > 6 GeV . 

To study the contribution from diffractive events, the 
PYTHIA 17 and PHOJET [H] Monte Carlo generators 
were used. Non-diffr active inelastic pp events were mixed 
with single and double diffractive events in accordance 
with the corresponding generator cross-sections for such 
processes. The diffractive contribution was found to be 
below 1% for the d(N)/dA<j> densities at p^ ad > 1 GeV in 
PYTHIA, and almost entirely concentrated at low multi- 
plicities (fewer than four topological clusters) . The contri- 
bution of diffractive events is larger for (d 2 N / drjdcf)) and 
(d 2 Pt /drjdcf)) measured at p!p d < 3 GeV , but becomes 
negligible for pif ad > 3 GeV . Diffractive contributions are 
higher for PHOJET, but their contribution was found to 
be smaller than the systematic uncertainties on the final 
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measurements. No attempt to subtract diffractive events 
from the final measurements was made. 



7 Systematic uncertainties 

The systematic uncertainties on the measured distribu- 
tions were determined by changing the selection or the 
analysis procedure and repeating the analysis. The largest 
uncertainties are described below: 

• The following procedure was used to estimate the ef- 
fect of the relative energy-scale uncertainty on the final 
measurements. The double ratio (E /p) mc / (E/p) was 
calculated for isolated single particles as described in 
Sect. [5] The effect of the energy-scale uncertainty on 
the measured densities was found by decreasing and 
increasing the pt of topological clusters in the Monte 
Carlo simulation, keeping the same cluster j?t in the 
data. The magnitude of the variation was set by the 
value of the double ratio calculated in a grid defined 
in r) and p. To simplify the calculation of the sys- 
tematic uncertainties, a common variation was used 
for all topological clusters independent of their ori- 
gin (hadronic or electromagnetic). The effect of the 
energy-scale uncertainty is significantly larger than that 
due to the event selection (including trigger) [15] . 

• The dependence of the bin-by-bin corrections on the 
detector material description was estimated by recal- 
culating the corrections using two further samples: one 
with an extra 10% of material in the tracking system, 
and the other with ~ 15% additional material in the 
region \ij\ > 2. 

• The physics-model dependence of the bin-by-bin cor- 
rections was estimated using the PerugiaO tune [22] 
instead of PYTHIA MC09. This uncertainty was sym- 
metrised. 

• A comparison of multiplicities of topological clusters in 
bins of track multiplicities indicated some discrepancy 
between data and Monte Carlo for events with low 
track multiplicities (see Sect. [3]). To estimate a system- 
atic uncertainty to account for this discrepancy, the 
bin-by-bin acceptance corrections were calculated af- 
ter re- weighting the PYTHIA MC09 detector-level dis- 
tributions. For this, cluster multiplicity distributions 
were measured in bins of track multiplicity and re- 
weighting factors were calculated by taking the ratio 
of the above distribution in data and PYTHIA MC09. 
The re-weighting procedure also addresses the uncer- 
tainties on the noise description used in the Monte 
Carlo simulation and other effects related to the cluster- 
reconstruction efficiencies. 

• A systematic uncertainty was estimated to account 
for differences in the pt resolution of leading topo- 
logical clusters in data compared to the Monte Carlo 
expectation. Discrepancies in the tails of the distribu- 
tions of p^ ad (dusters) /p^ ad (tracks) were used to ex- 
tract weighting factors, which were then used to recal- 
culate the acceptance corrections. 



Table [T] shows the values of the systematic uncertain- 
ties discussed above as a percentage of the measured val- 
ues. Only the largest values are shown for the bins with 
the most significant effect from the selection variations or 
change in the experimental procedure. 

In addition to these uncertainties, the following sys- 
tematic variations were also included: 1) in order to re- 
duce the contribution from diffractive events, the mea- 
surement was repeated after removing events with fewer 
than four clusters; 2) the positions of cluster centres in rj 
and (f> were shifted by the size of one cell; 3) an alternative 
model (FTFP-Bertini) for the hadronic-shower simulation 
in Geant4 was used to extract the correction factors; 4) 
the calorimeter transition region of 0.94 < |r/| < 1.06, 
which is not well described by the Monte Carlo simulation 
was removed in the data and in the simulated PYTHIA 
MC09 sample. These variations each give systematic un- 
certainties below 2%, with the exception of that for diffrac- 
tive events which indicate a 7 — 10% systematic uncer- 
tainty for the (d 2 N / dr]d(j)) and (d 2 Pt / dr]d<p) densities 
measured at p^ ad < 3 GeV . As an additional systematic 
check, the measurement was also repeated using topolog- 
ical clusters at the electromagnetic energy scale and simi- 
lar differences between data and Monte Carlo simulations 
were observed. 

The overall systematic uncertainty was determined by 
adding the above uncertainties in quadrature. 

8 Results 

Figure[5]shows the density distribution d(N) / dA<fi of stable- 
particles as a function of the distance in azimuthal an- 
gle between the leading particle and other particles in an 
event for ^/s = 900 GeV . This density, defined in Sect. El is 
calculated for events selected by requiring a particle with 
p}f ad above the values indicated on the figure. The detec- 
tor correction for this density is discussed in Sect. El The 
total uncertainty, computed from the addition of statisti- 
cal and systematic uncertainties in quadrature, is shown 
as a shaded band on all measurements. 

The angular distribution shown in Fig. [S] has a peak 
at A<fi — which reflects the particle activity from the 
hard interaction. The peak narrows as plf. ad increases. The 
shape of this distribution is similar to that observed in 
the recent track-based publications , and also similar 
to the transverse-momentum flow around jets observed 
at a lower pp collision energy |34) . The particle densities 
measured using topological clusters are higher than the 
charged-particle densities measured using tracks [5] , which 
is expected from the neutral-particle contribution. 

Figure IHl shows the A(f> density distributions for y^i = 
7 TeV . The distributions show narrower peaks, for a given 
p^ ad threshold, than for the y/s = 900 GeV data. 

The data are compared to the PYTHIA Monte Carlo 
generator with the MC09, PerugiaO and DW tunes, PHO- 
JET and HERWIG+ JIMMY. The Monte Carlo genera- 
tors reproduce the general features of the data, but fail to 
describe the detailed behaviour, as can be seen in the fig- 
ures. The MC09 and PerugiaO PYTHIA tunes are closest 
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Check 


d{N)/dA<t> 


(<fN/drid4>) 




Energy scale 


±4.3% 


±4% 


±5.6% 


Additional material 


+3.5% 


±3% 


±3.6% 


Model dependence 


±3.5% 


±5% 


±4.5% 


Multiplicity reweighting 


±4.5% 


±10% 


±11% 


Resolution reweighting 


±0.4% 


±6% 


±6% 



Table 1. A summary of the most important systematic uncertainties. The table lists the values of contributions from different 
groups of systematic checks. Only the largest values are shown, taken from the bins with the largest effect when the systematic 
variation was applied. 




Fig. 8. The average number of particles per unit of pseudorapidity as a function of the azimuthal separation between the leading 
particle and other particles in inelastic pp collisions at y/s = 900 GeV . The shaded band shows the statistical and systematic 
uncertainties added in quadrature. The densities are shown for (a) ptf ad > 1 GeV , (b) p£ ad > 2 GeV and (c) p^ d > 3 GeV . 




Fig. 9. The average number of particles per unit of pseudorapidity as a function of the azimuthal separation between the leading 
particle and other particles in inelastic pp collisions at y/s — 7 TeV . The shaded band shows the statistical and systematic 
uncertainties added in quadrature. The densities are shown for (a) p l ^ d > 1 GeV , (b) p^ d > 2 GeV and (c) p^ d > 3 GeV . 
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Fig. 10. The average number of stable particles per event per unit interval in T)-(j>, as a function of p^ ad , for the transverse 
region indicated in Fig. [TJ The results are shown for (a) y/s = 900 GeV and (b) for y/s — 7 TeV pp collisions. The shaded band 
shows the statistical and systematic uncertainties added in quadrature. 



to the data. The PHOJET generator significantly under- 
estimates the particle densities, while the PYTHIA DW 
and HERWIG overestimate the data at A(f> ~0. The data 
are seen to have a large discriminating power and are thus 
useful to constrain the parameters of Monte Carlo gener- 
ators. 

Figure [TU] shows the mean number of particles per 
event per unit interval in r\ and <p as defined in Sect. [6] The 
density was calculated in the transverse region illustrated 
in Fig. 12 as a function of p!p d . None of the Monte Carlo 
predictions describe the data well. The DW tune is the 
most similar to the observed data. As is seen in the A<j> 
distribution, the PHOJET simulation lacks a hard com- 
ponent for y/s = 7 TeV . The particle density increases 
almost by a factor of two, going from y/s = 900 GeV to 
y/s = 7 TeV at a similar p!f ad , which is comparable to what 
is seen in all Monte Carlo generators. 

Figure [11] shows the mean scalar p^ sum for stable 
particles in the transverse region as a function of pl£ ad . As 
for the particle densities, the mean transverse-momentum 
sum is measured per unit interval in 77 and (f> (see Sect. [5]). 
Again, the Monte Carlo predictions do not fully describe 
the data. The largest discrepancy with the data is found 
for the PHOJET generator. 



9 Conclusions 

Particle densities sensitive to the underlying event in pp 
collisions at centre-of-mass energies of 900 GeV and 7 TeV 
are presented. This is the first such analysis completely 
based on calorimetric measurement of three-dimensional 
energy depositions, which is made possible by the fine 
granularity of the ATLAS calorimeter with transverse and 
longitudinal samplings. 

The particle densities were studied and compared with 
several Monte Carlo generators tuned to pre-LHC data. 
None of the Monte Carlo generators describe the mea- 
surements well. In particular, the Monte Carlo predictions 
have discrepancies with the data for the particle density as 
a function of the azimuthal angle between the leading par- 
ticle and any other particle in an event. The Monte Carlo 
generators typically predict a lower particle density in the 
transverse region (\A<p\ ~ t/2), while in the toward re- 
gion (Acf) ~ 0), the PYTHIA DW and HERWIG+JIMMY 
generators both overestimate the densities. PHOJET sig- 
nificantly fails for the y/s — 7 TeV data. For the particle 
densities as a function of p^ ad , all the Monte Carlo gener- 
ators also fail to describe the data, predicting lower than 
observed particle activity in the transverse region. A sim- 
ilar conclusion holds for the total transverse momentum 
of particles in the transverse region. 

The particle densities measured using topological clus- 
ters are higher than the charged-particle densities mea- 
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Fig. 11. The average scalar pt sum for stable particles per unit area in rj-(j) in the transverse region as a function of pl£ ad 
for (a) yfs = 900 GeV and (b) y^s = 7 TeV data. The shaded band shows the statistical and systematic uncertainties added in 
quadrature. 



sured using tracks |5j. This is expected from the neutral- 
particle contribution. The discrepancies between the data 
and Monte Carlo generators agree with those observed for 
charged particles [5 . These measurements have system- 
atic uncertainties independent to the track-based studies 
and provide additional information which may be used to 
improve the Monte Carlo description of the complete final 
state produced in pp collisions. 
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